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Abstract 

We describe how the matrix integral of Imbimbo and Mukhi arises from a limit of the 
FZZT partition function in the double-scaled c = 1 matrix model. We show a similar 
result for OA and comment on subtleties in OB. 



1 Introduction 


Non-critical string theories with c < 1 are an important conceptual laboratory for exploring 
the inner workings of string theory. As world sheet theories, they are formulated by coupling 
c < 1 matter to Liouville theory and describe strings propagating in a space-time with a 
linear dilaton and a tachyon condensate. While these world sheet theories are, in general, as 
difficult to solve as critical string theories, it is a special feature of low-dimensional strings 
that there often exists a dual description in terms of an exactly solvable matrix model [1-5]. 

There are, typically, two distinct formulations of a given c < 1 string as a matrix model: 
the double-scaled matrix model and the Kontsevich model. Historically, these models were 
derived in completely different ways. The double-scaled matrix model is based on ‘t Hooft’s 
idea that, at large N, the Feynman diagrams of gauge theory can be thought of as random 
triangulations of a string world sheet. The Kontsevich model, meanwhile, is found by directly 
solving the flow equations implied by the integrable structure. Amazingly, this solution takes 
the form of a simple matrix integral [6]. Both approaches can be used to compute observables, 
such as correlation functions of vertex operators on the world sheet, and give identical results. 
However, until recently, it was not understood how they are related to each other. 

The existence of dual formulations, one in terms of a string theory and the other in 
terms of matrices, is reminiscent of open/closed-string duality and suggests that the matrix 
models may arise from open-string theories living on branes. There are two classes of branes 
in low dimensional string theories, which arise from two possible boundary conditions in the 
Liouville sector: FZZT branes [7,8], which are the analogue of Neumann boundary conditions 
and are extended, and ZZ branes [9], which are the analogue of Dirichlet boundary conditions 
and are localized in the strong coupling region. 

ZZ branes are unstable. This, together with the fact that they live in the strongly 
coupled region, makes a direct construction of the open string theory living on them difficult. 
However, it was shown in [10,11] that the decay of a ZZ brane can be interpreted as a single 
eigenvalue of the double-scaled matrix model rolling down the unstable potential. This allows 
one to think of the double-scaled matrix model as the open string theory living on an inhnite 
number of ZZ-branes that have condensed to form a Fermi-surface. 

FZZT branes are stable. They are labeled by a parameter /i^, known as the boundary 
cosmological constant. Although the FZZT brane is the analogue of Neumann boundary 
conditions, it is only semi-inhnite in extent. For example, if we denote by \hb) the boundary 
state of an FZZT-brane in the {p, q) minimal model, then, in the mini-superspace approxi- 
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0 = (1/6) log(/iB) 




Figure 1: Illustration of profile ~ e of the boundary state \^b) 

mation, we have the prohle^, 

(0|/iB) ~ where b^ = -, c=l- . (1.1) 

q pq 

As one can see in hgure d the brane exists throughout the weak coupling region at large 0, 
but disappears around 0 ~ -^logps. Hence, by varying ps, FZZT branes can be used as a 
probe of target-space [12], 

Having identihed the open string theory on ZZ branes with the double-scaled matrix 
model, it is natural to try to identify the Kontsevich model with open strings living on FZZT 
branes. The precise relation was hrst pointed out by Gaiotto and Rastelli who considered 
Witten’s cubic open string held theory [13] for the open strings ending on stack of n FZZT 
branes for the case of the minimal {p,q) = (2,1) model [14]. Remarkably, they found that 
the open string held theory [13] on a stack of FZZT branes could be reduced to precisely the 
action of the n x n Kontsevich matrix integral^, 

R(M, A) = TrMA-. (1.2) 
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Unfortunately, while this result is elegant, it has not been easy to generalize it to other 
examples, such as the {p,q) ^ (2,1) models [14,15]. 

An alternate approach, which seems to generalize more easily, is to compute the open 
string theory on the FZZT branes using the double-scaled matrix model. An FZZT brane 
can be incorporated into the double-scaled matrix model by inserting into the path integral 
the exponential of the loop operator, 

exp(Trlog(/iij - $)) = - 4*), (1.3) 

where $ is the N x N matrix of the double-scaled matrix model. This approach was pursued 
in [16], and it was shown that one could derive the Kontsevich model from the partition 
function of the (2,1) matrix model in the presence of a stack FZZT branes. 

^Our convention is that (j) +oo is the weakly coupled region. 

^Sometimes, the action is presented in a form where the term linear in M is replaced by term quadratic 
in M by shifting M. 
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Unlike the string field theory analysis, the donble-scaled matrix model analysis can be 
generalized in a straightforward fashion to other examples. For example, it was shown in [17] 
that, using the the double-scaled {p, 1) model [18] in the presence of n FZZT-branes, one 
could derive the corresponding generalization of the Kontsevich integral. 


Z(A) 


dMe 



MPjp_\ 

P+1 J 


(1.4) 


In spite of the success of these arguments in c < 1 models, the extension to c = 1 has 
not been worked out. As with the c < 1 models there is an analogue of Kontsevich matrix 
integral, 

Zn{A,t) = {detAY J u = -ip , (1.5) 

first derived by Imbimbo and Mukhi [19]. However, the connection between the open string 
theory on FZZT-branes and the Imbimbo-Mukhi model has not been made. 

Several attempts to make this connection can be found in the literature. One approach is 
to work out the string field theory living on a stack of FZZT-branes in the same spirit as the 
Gaiotto-Rastelli computation. It is not clear, however, how the degrees of freedom of open 
strings ending on a semi-infinite world volume should arrange themselves into parameters A 
and t of the Imbimbo-Mukhi model and the authors of [20] did not succeed in recovering the 
Imbimbo-Mukhi model this way. 

Another method of deriving the Imbimbo-Mukhi model used the so-called normal-matrix 
model [21]. This model is similar to the Imbimbo-Mukhi model in that it is a matrix 
integral solution to the integrable flow equations. As discussed in [22], one can perform 
a computation similar to [16,17] in which one evaluates the expectation value of certain 
determinant operators in the model to reproduce the Imbimbo-Mukhi model. While this 
gives a nice connection between the normal-matrix model and the Imbimbo-Mukhi model, 
the relationship between the determinant operators in the normal-matrix model and FZZT 
branes is not clear. 


2 FZZT branes in bosonic c = 1 string theory 

Having reviewed what is known about the c = 1 case, we can now turn to the main goal of 
this paper: we wish to fill in the conceptual gap between the open string theory on the FZZT 
branes and the Imbimbo-Mukhi model. To do so, we attempt to compute the expectation 
value of the exponential of the loop operator, as was done in the c < 1 case, and manipulate 
it into the form of the Imbimbo-Mukhi matrix integral [17]. As we will see, many of the 
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ingredients that were used in the c < 1 computation have analogues in the c = 1 case. There 
is, however, one new subtlety. As described below, one must take a certain scaling limit of 
the FZZT partition function before one can derive the Imbimbo-Mukhi model. 


We begin our discussion with a brief review of the relevant aspects of the c = 1 matrix 
model that we need for our derivation. We mostly follow the approach of Kostov [23,24]. 
The matrix-model action is given by 


where 


S 



* (-i(V^4)''+ p>"), 




( 2 . 1 ) 

( 2 . 2 ) 


and <h(a:) and A{x) are N x N Hermitian matrices. For the bosonic theory, we work in units 
where a' = 1. Since we are interested in deriving the Imbimbo-Mukhi model, we work with 
the Euclidean action and compactify the time direction with period [3 = 2ttR. 


It is useful to rewrite the action, (EH), in chiral form. To do this, one hrst switches to 
the hrst order form, 

rP 

S = Tt dx (*P(VA<h) - \P^ + , (2.3) 

Jo 

which reduces to equation EH upon integrating out P. Then, substituting 

$±P 


= 




one arrives at the chiral action. 


rP 

S = Tt dx{-iX+VAX. + X+X_). 
Jo 


(2.4) 


(2.5) 


In this form it is relatively straightforward to reduce the path integral over X±. and A to an 
ordinary three-matrix integral [23-25]. Dehning X± = X±{t = 0) one hnds. 


Zjv = / dX+dX.dne‘'^f’‘*’‘—’ ‘), 


( 2 . 6 ) 


where the subscript, N, is a reminder that the partition function depends on the rank, N, 
of the matrices and q = e*^. The matrix, hi G U{N), is given by 

dio MAdx 


n = p 


(2.7) 


This expression can be simplihed further by integrating out hi, which reduces the partition 
function to an integral over the eigenvalues of X± , 


N 


Zn — 




k=l 


) det(e 

jk 


iq ^X. x+ 


( 2 . 8 ) 
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In order to extract the observables of the c = 1 matrix model, one compntes the grand 
partition function, which, following [24,25], can be expressed in the form of a Fredholm 
determinant. 


Z{fi) = = det(l + e-'’>‘K_K+), 


(2.9) 


where 

lK-f](x+)^ fdx_e-'^*’-f(x_) ( 2 , 10 ) 


are operators acting on functions of and X-, respectively. The Fredholm determinant 
is computed by diagonalizing the operators and K_ and taking the product of their 
eigenvalues. This closely resembles what one does in computing the partition for the two- 
matrix models in terms of orthogonal polynomials. Up to non-perturbative corrections, the 
partition function of the undeformed Euclidean c = 1 model can be written as 


logZ(/i)= i(j)o{-n + 

r=Z+| 



^^ME) ^ ^iW4 Cosh(7rU) ^ 
TT 



( 2 . 11 ) 


As is well known, when compactihed on a circle, the c = 1 theory contains a tower of 
tachyon deformations with discrete Euclidean momenta. These deformations change the 
asymptotic shape of the Fermi surface in phase space. One can incorporate them into the 
matrix model by generalizing (12.8p to 

Z„ = / d<|t(e-V<) d|t(e‘''"V4), (2,12) 

^ fc=l 


where 

[dx±\ = U±{x±) = iz±kX^J.^ (2-13) 


and the Zk are the deformation parameters corresponding to tachyons with momentum k. 
The sign and the phase of Zk was chosen so that it coincides with the deformation parameter 
of [26]. This identihcation is explained in more detail in appendix A. The grand partition 
function is still given by the Fredholm determinant, but for the deformed operator. 


Z{ii) = det(l + e'’'‘/f_/f+), 


(2,14) 


where 

[fun(x-) = j [*+]e‘’"‘"+-/(i^+). 

lK-f](x+) = j [da:_]e-“+'^-/(a:-) . (2.15) 
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It was shown in [23] that snch a deformation can be described in terms of a Toda flow and 
that the partition function is related to the r-function via^ 

Z(p) = r - |) + |,z+,z_) . (2.16) 

The initial condition for the flow, (ITTTl) . can also be expressed in the form of the string 
equation, found in [30,31]. 

Having given a brief review of the c = 1 matrix model, we now study what happens 
when we add an FZZT-brane to the story. As mentioned earlier, an FZZT-brane can be 
incorporated in the matrix model by inserting into the path integral an exponential of the 
macroscopic-loop operator, 

= det(<|)(x) - = det^(X+(a:)X_(a;))/A/2 -/isj , (2.17) 

as hrst considered in [32]. Note that the determinant operator generates an FZZT brane with 
Dirichlet boundary conditions in the c = 1 direction and, hence, is localized in Euclidean time 
at a point x. The “±” in the exponent is included so that we can simultaneously consider 
the determinant and the inverse determinant operator insertions. As we will see later, the 
connection with the Imbimbo-Mukhi model favors the inverse determinant [22]. However, 
inserting the determinant also gives rise to a matrix integral that encodes the r-function 
of the Toda integrable flow and is, roughly, the complex conjugate of the Imbimbo-Mukhi 
model. 

Inserting (ITT7I) into the chiral Lagrangian, (ESD, one sees that the determinant has a 
simple X dependence, 

^±w{x) ^ det(^(e“/^A:+(t = 0) + e-“/^A:_(t = 0))/v^ - (2.18) 

One can also see this x dependence by considering the Heisenberg form of the operators, 
X±, in the Hamiltonian language. 

Now, it would be nice if inserting this operator into the path integral was equivalent to a 
tachyon deformation. It clearly is not, since the determinant contains a complicated mixture 
of A+ and X_ operators, and does not factorize into an term and an X_ term in the 
way that the tachyon deformation, (EH, does. However, as we now show, it is possible to 
get a pure tachyon deformation if we take a certain limit of the determinant operator. 

Suppose that we analytically continue x = —it and send t > cxd while simultaneously 
rescaling 

fiB = eVn • (2.19) 

^The connection between c = 1 matrix model and the Toda integrable hierarchy was first pointed out 
in [27]. For introduction to Toda lattice hierarchy, see [26,28,29] 
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Then, up to an overall normalization factor, the operator asymptotes to 

det -. (2.20) 

In the case where the radius of the c = 1 direction is one, inserting this operator into the 
partition function is equivalent to a shift in the Toda flow, 

5zk = 5k > 0 . ( 2 . 21 ) 

n 

Readers familiar with [22] will recognize this deformation as the same one generated by one 
of the determinants in the normal-matrix model, det(l — Z/ The other determinant, 
det(l — Zy I can generated by taking the opposite sign in the analytic continuation, 
X ^ it . Notice that away from R = 1 there does not seem to be any relation between the 
FZZT-brane determinant and the determinant of the normal matrix model. 

Once we have identified that our limit of the FZZT determinant generates the Toda flow, 
the integrable structure of the model is enough to determine that the expectation 
value of will give the Imbimbo-Mukhi model, as was shown in [22] using the normal- 

matrix model. However, we would like to understand how to derive the Imbimbo-Mukhi 
model more directly. 

As in the c < 1 case, it is useful to introduce the Baker-Akhiezer function, T^(A), which 
is given in terms of the r-function by 

^ llf+ZikATizAexp[z+(A'''*) + <ilog(A‘''«)], 

T\Zj^^ 2 '_, S) 

>I-;(A) = h£±li_ZiAAl£±Fexp|z_(A-‘/«) + s log(A‘'">)], (2.22) 

where 

/ a2 \3 \ °° 

e(A)= ( A, y,, z+(A) = ^x„A", x_(A-') = ^ g_„A-" . (2.23) 

For a formal definition of the Baker-Akhiezer function in terms of dressing operators in the 
Lax formalism, see equations (25) and (27) of [33]. 

Notice that, in the expression for T+, the argument of the r-function is shifted by 
2 + — e(A“^/'^). Setting A = /ig and i? = 1, we see that this is the same shift as (1^ . 
provided we consider the inverse determinant. Hence, the r-function in the presence of the 
inverse determinant, (ESOD, can be computed if we know the Baker-Akhiezer function. 

By itself, this fact is not very useful. However, the Baker-Akhiezer functions are also 
closely related to the biorthogonal eigenfunctions, of fhc deformed operators, K± given 
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in ()2.15|1 [23]: 




s=-R{iE+l/2) 




s=-R{iE+l/2) 


(2.24) 


Since the biorthogonal eigenfunctions are related to each other by Fourier transform, it 
follows that and 4/7 are related by 





(2.25) 


Now, suppose that = 0. Then the non-trivial pieces of drop out and is completely 
determined as a function of . 2 _ and A; 


T, (A_) = A/'exp[z_(A_^^-^) + slog(Ai^^)], 


(2.26) 


where A/” is a normalization factor independent of A and z^, which we can ignore. Changing 
variables to m = Al^, one has 


=jdn. exp|«m + - (s + K) log(mVH)|. 

Using (ESI, one hnds 

I dm (mar'-eMiam + ._(»■/«) _ 

which is equivalent to 

Zi(A, t_) = r(-e(a-^/^), Z-,s) = J dM (2.29) 

with the identideations, 

t_n = z^n, fn = a = s = i^R — —, u = —ifi. 

(2.30) 

As we discuss in appendix C, this argument can be iterated to give the general expression. 


Zn{A,tJ) = / dMi 


A{M) 




ixi?+(i/-l)/2g-ix((MiAi)^+U-nMf)-nlog(Mi) 


with Miwa-Kontsevich-like transform. 


t 


n 



(2.31) 


(2,32) 


This is the same result obtained in [22] using the normal-matrix model.Setting i? = 1, we 
reproduce precisely the Imbimbo-Mukhi integral. 

^It should be emphasized that while the sign of z„’s is fixed by matching to the conventions of [26], there 
is no intrinsic definition for the sign of the t„’s. In particular, the sign convention for the t„’s in [19] and [22] 
appear to differ by a sign. 








Recall that, to match with the argument of the Baker-Akhiezer function, we had to 
consider the inverse determinant. This is in agreement with the observation made in [22] 
that the Imbimbo-Mukhi integral is derived from the insertion of an inverse determinant in 
the normal-matrix model. 


One can also consider the insertion of a determinant. In this case, one must work with 
the dual Baker-Akhiezer functions [29]^, 


*r(A) = 


r(z+ + e(A z_,s) 


exp 




’■(z+.z-.s) 

As with the Baker-Akhiezer functions, they are related by Fourier transform; 


1 


from which one can derive. 


T+*(A+) = / dA:^e-*^+^-'^Ay^T;"(A_), 




(2.33) 


(2.34) 






—_ A-k 


V = Ifl 


(2.35) 


For i? = 1, one can use the Harish-Chandra/Itzykson-Zuber formula [34,35] to re-express 
this integral in matrix form. This again yields a matrix model identical to the Imbimbo- 
Mukhi model. The only difference between and (lOTl) is that we have changed the 

relative sign between v and z/i. Since both models are an even function of /i, this does not 
have any effect on the closed string observables. Therefore, at the level of reproducing the 
r-function, one can not distinguish between the determinant and the inverse determinant. 
This choice is ultimately tied to whether the FZZT-ZZ strings are bosons or fermions. 


When i? ^ 1 it is still true that (E3T1) and give complete information about the 

r-function. However, since the FZZT determinants no longer give rise to a Toda flow, (jZSH) 
and (IT3IK11 are no longer related - in any obvious way - to FZZT-branes. Thus, we hnd that 
the open/closed string duality, relating the Kontsevich-like matrix integral and the FZZT 
brane, appears to work only at self dual radius. This is consistent with the observation 
in [22] that (j2.3Hl and (j2.35jl cannot be represented in the form of a matrix integral except 
at R — 1. After all, if fj2 .31 jl and were derivable from an open/closed duality, one 

would expect the result to take the form of a matrix integral for general R. 

^They are related to “ )*(s, z+,Z-,X) in the notation of [26] according to 
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3 FZZT branes in c = 1 string theories 


Our analysis relating the determinant operator to a Kontsevich-like matrix integral can be 
extended to the c = 1 models [36-38]. There are two c = 1 models: type OA and type OB. 
We consider them separately. 

3.1 Type OA 

The type OA theory is formulated in terms of a t/(A) x U{N+q) gauged matrix model [37,38]. 
The parameter q gives the number of units of Ramond-Ramond flux in the background.® As 
in the bosonic case, one can fix the gauge and reduce the matrix degrees of freedom to the 
eigenvalues. This gives rise to a model which is equivalent to the deformed matrix model 
originally constructed in [40]. We will work in units where a' = 1/2 so that the harmonic 
oscillator has the same oscillator frequency as the c = 1 model. Physical observables such 
as the S-matrix and the partition function of this model were computed in [41]. These 
computations were then reformulated in light-cone variables [39], which are better suited for 
the analyzing the deformations and Toda integrable structure of the model [42].^ 

The main difference between the type OA model and the bosonic c = 1 model is that, in 
OA, the wave functions, and 'ip-ixJ), are related by a Bessel transform instead of a 

Fourier transform. As a result, the Baker-Akhiezer functions are also related by an integral 
Bessel transform, 

*yA+) = r«:' \/a+a:V,(a+a;‘)a_®;(a_), (3.1) 

The partition function can be written in the form®, 

„ N 

Zn= W[dxl][dx^] det Ux]x^Jq{xxxl) \ det UxjxlJx]xl)\ , (3.2) 

d k=i ^ ^ 

where 

[dx±\ = U±{x±) = iz±kX^J.^. (3.3) 

We can also express Z{^) as a Fredholm determinant, 

Z(/i) = det(l + e^^J?_A+), (3.4) 

®See [39] for a discussion of the subtleties associated with the two ways of introducing charges in this 
model. 

^See also [43-46] for a related discussion in the context of minimal models. 

®Note that q here refers to the background flux and not e®^. 
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where 


lK*f]{x-) 

lK-f](x+) 


j [dx^]^x+x-j,(e ''’x+x-)f(x+) , 
j[dx.]^x+x_J,(-x+x-)f(x-) . 


(3.5) 


Using the relation between the r-fnnction and the Baker-Akhiezer fnnction, , one 

hnds 

r(=Fe(a“^'^^), z_, s) = J Jg{am) exp[±z_{m^^^) — \og{m)], (3.6) 

which generalizes (IT^ to type OA. Just as in the c = 1 case, one can rescale 

= m = a = ii^l'^A^ (3.7) 

to bring this expression into an Imbimbo-Mukhi form. Generalizing to the multiple FZZT 
case one hnds the expression, 


r(=Fe(a ^/^),z_,s) 

A(M) 


= / dM 


n(M,A)-+^/V,(±/iAfMf) exp[±*p(t_(M,)) - nlog(M,)], (3.8) 


where 


Using 


s = ±f/ii? — 


R-\ 




—uma cos 9 


= J,(m) , 


(3.9) 


(3.10) 


one can integrate in an additional variable and express the shifted r-function in a two-integral 
form; 


r(=Fe(a ^/^),z_,s) = 


n 


v/2+(R-l/2)/m+q/2 


POO ^ / -j 

I 


-SiZi+q 


A(z)i/2« 11 

/ n 

<i!/,A(s)‘'“ n!/r'"+'''-‘-'">'“-<"+«'''>£!-«/‘-exp|->.((_(!/.'''"®))] , (3.11) 

i=l 

where we have used the change of variables. 


m 


Hi = rrij , s = — cos 0 , Zi = , u = . 

au 

Setting R = 1/2, one recovers a matrix integral representation found in [47,48]. 


(3.12) 
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To relate the Baker-Akhiezer function and the r-function to FZZT branes, we consider 
the determinant operator [38], 

^±w(x) ^ det(X(a;)X(x)^ — = det ^(X+(x) + X_(a;))(X|(x) + x !_( x ))/2 — jj,^ 

(3.13) 

Taking the same scaling limit as we did in the bosonic case, (nmn) reduces to 

det ( 1 - j = . (3,14) 

Comparing this with the parameterization of the deformation parameter, (ini, we see that 
the exponentiated macroscopic-loop operator matches with the Baker-Akhiezer function at 
R = 1/2. This provides the link between FZZT branes and the Imbimbo-Mukhi integral, 
(EH), as well as the matrix integral of [47,48]. 

Since c = 1 models are on firmer non-perturbative footing that the bosonic c = 1 models, 
it is instructive to examine the the deformed r-function function at the non-perturbative 
level. To that end, it is useful to go back to the integral transform, (EH), and write it in the 


form. 

/a+a:* (ff+(A+A:') + h-{x^\-J)) >i<1(a_), 

(3.15) 

where 


J^±(^) _ Jq{x+X-)±iYg{x+X-) 

(3.16) 


are the Henkel functions. Assuming that the background is only deformed by .s_, one can 
use the relation between the Baker-Akhiezer function and the r-function to write 


r(—e(a Z-, s) = / dm{ma) — 

Jo 

(3.17) 

To dehne the theory non-perturbatively, one must specify the contour of integration for the 
variable, m. For the undeformed background, = 0, one can integrate the Henkel function 
H^{x) anywhere along the upper/lower half of the complex x plane. This choice of contour 
reproduces the correct perturbative expansions in 2 ;_. 

When the deformation parameters are finite, one encounters a discrete choice of al¬ 
lowed contours that give rise to finite but distinct results. To illustrate this point, consider 
the simplest case of i? = 1/2 with z^i being the only non-vanishing deformation parameter. 
We will also take p to be real and positive, and t_i to be real but have small negative imag¬ 
inary part. Then, z is mostly imaginary with a small negative real part. The asymptotic 
behavior of the integrand (IXT7I) is dominated by the exp{z-imJ) factor and converges when 
integrated along contours that go to inhnity along the unshaded region, as illustrated in 
hgure m 
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A 

Re[m] 


Figure 2: An illustration of possible contours of integration in m-space for which the integral 
(IXT7I) converges. Contours asymptoting to infinity along the shaded regions are excluded. 
Here it is assumed that z^i is the only non-vanishing deformation parameter. 

The integral involving H^{am) can be performed along either of the contours labeled 
by A and B in the hgure. The difference between the two contours will only contribute 
non-perturbatively in the asymptotic expansion in small For a more 

general deformation, there could be many distinct contours consistent with the perturbative 
expansion. Therefore, in order to dehne the r function non-perturbatively, we need a definite 
prescription for the contour. 

One natural requirement is that the contour have a well dehned limit as the 0. 

Let us assume, for simplicity, that there is some k > 0 such that t_n = 0 for n > k.^ Then, 
the factor of exp[z_(m^/^)] is always oscillating along the positive real axis. By adding a 
small imaginary part with appropriate sign to one can open a small region around the 
positive axis for which the exponential factor exp[z_(m^/'^)] converges. This is illustrated in 
figure 121 Both Henkel functions can then be integrated over contours on the correct side of - 
and arbitrarily close to - the real axis. In figure |21 this corresponds to choosing the contour 
A. We propose that this contour is the correct non-perturbative dehnition of our integral, 
(uni). This prescription gives a natural, though not unique, non-perturbative dehnition 
of the Toda how. It would be instructive if one could hnd some consistency check for this 
prescription. 

®Such deformations were also discussed in [49]. 
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3.2 Type OB 


Finally, let us briefly describe the extension to the case of type OB theory. This model turns 
out to be more subtle at the non-perturbative level than the type OA model. 

We hrst consider the case of vanishing background RR flux. The model is formulated 
as the quantum mechanics of an inverted harmonic oscillator potential with fermions hlling 
both sides of the potential to the same level. There is an explicit symmetry exchanging 
left and the right side of the potential, and, as long as the unbounded potential is regularized 
in a way that preserves the left/right symmetry, the spectrum of states can be separated into 
sectors which are odd and even with respect to this Z 2 . Since these two sectors are decoupled 
from one another, the partition function factorizes into a product of two r-functions [42], 


•^(h) '^oddi^^odd+1 ^odd—1 ^'}T~even(^^even+: 




(3,18) 


The Baker-Akhiezer functions for the odd and even sectors have an integral transform, 

poo 

^W(A+) = / dX-J sin(A+A:^)vl/:,,_(A_), 

poo 

«'™„+(A+) = / d\-J cos(A+A;')4')„„_(A_). (3.19) 

Jo 


Using the relation between r-functions and the Baker-Akhiezer functions, one arrives at 
an integral expression of the form, 

fc) 

A(M) 


dM- 


cosh(,.AlfM«) exp|rf_(M,) - n log(Mi)], 


^oddi^t'k') t—k') '^) 


dM, sinh(i/AfMf) exp[ut_{Mi) - nlog(M,)] , (3.20) 

^ ' 1=1 


with 


tk 


A 


-k 


uk 


(3.21) 


These expressions are to be thought of as arising from the Fredholm determinants, 

Teven = det(l 6^^ K^^en-Keven+), Todd = det(l 6^^ Kodd-Kodd+), (3.22) 


where, for example. 


[^even+ 

j"[dx+]cos{e *^x+x_)/(x+). 


[Kodd+f]{x+) = 

1 [dx_] sm{e-^^x+x_)f{x_), 

(3.23) 
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with [dx±\ defined in the same way as in ()2.1()|1 . 

To relate these expressions to FZZT branes, observe that the OB FZZT determinant 
[36,44,50], 

g±w(x) ^ ciet(X(x)^ — = det ((X+(a:) + X_{x)Y — , (3.24) 

asymptotes to 

det (^1 - , (3.25) 

which, for R = 1/2, matches the deformation of the tk implied by the Baker-Akhiezer 
function. In a broad sense, we have therefore succeeded in relating a scaling limit of FZZT 
amplitude to a Kontsevich-like integral expression dsm. Note however that unlike in the 
bosonic and the OA cases, R = 1/2 is not the special radius for which (I3.20j) can be re¬ 
expressed in matrix form using the Harish-Chandra/Itzykson-Zuber formula. 

Note, also, that the exponentiated macroscopic-loop operator perturbs both the even and 
the odd sectors simultaneously. As such, it does not probe the two sectors independently. A 
different object must be introduced to generate purely even or odd deformations. It is not 
clear what that object is. 

To make sense of the integral expression, ( 11201 ), at the non-perturbative level, it might 
seem that all one has to do is address the issue of contours as we did for the OA. However, 
there is an additional subtlety in the OB case. 

To understand the origin of this subtlety, consider turning on a flux. There are two 
proposals for the theory deformed by a flux. One method is to represent the flux by a 
relative shift between the chemical potentials of the even and odd sectors. However, it was 
found that such a prescription gives rise to violation of T-duality at a non-perturbative 
level [39,51]. 

An alternative prescription, where one adjusts the chemical potential of the left and 
right-moving sectors independently, was proposed in [39]. This prescription leads to partition 
function that is T-dual at the non-perturbative level. A second advantage of this formalism 
is that there is a well defined notion of a fermi-sea, which is lost when one deforms the 
even/odd sectors separately. 

Unfortunately, the integrable structure of OB has only been understood in the case where 
we consider a flow which acts only on the even or odd sectors separately. In this formalism 
it is possible to turn on a flux using the first method, in which the even and odd sectors 
have different chemical potentials, but it seems very difficult to consider the case where the 
left and right sectors are modified. Moreover, by splitting the fock space into even and 
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odd sectors, one loses the notion of a single Fermi-sea once any deformations are turned on. 
This suggests that Toda flow given in [42] may not represent a physical deformation of the 
background. 

Perturbatively, none of this is an issue since the difference between the left/right and 
even/odd bases only enters non-perturbatively. However, to make the connection between 
FZZT-branes and Konsevich-like integral expressions at the non-perturbative level, it seems 
one must hrst gain a better understanding of the integrable structure of 0B.^° 

4 Conclusions 

In this paper, we found a precise relationship between the Imbimbo-Mukhi model and the 
open string physics living on FZZT branes. The key ingredient in this relationship, which did 
not play a role in the c < 1 theories, was our scaling limit. Allowing the time coordinate of 
the FZZT brane to be complex, we pushed the location of the brane into the far Minkowski 
past. This reduced the complicated action of the full FZZT determinant to a pure Toda 
flow. 

Using analytic continuation to reduce a brane amplitude to a purely closed string process 
is reminiscent of the imaginary brane analysis of [52]. There one does not require a limiting 
procedure; rotating the brane to imaginary time is all that is required to decouple the open 
string states. This may well be true in our story; however, the closed string states would, in 
general, be a mixture of incoming and outgoing states. The extra step of pushing the brane 
to the far past is probably required to reduce the FZZT brane to a purely incoming source. 

We also considered the generalization to c = 1 models. For type OA theory, we found 
that the integral expression generalizing the Imbimbo-Mukhi formula [47,48] can be related 
to the scaling limit of a macroscopic-loop operator provided the radius is R = {l/2)\/2a'. 
We also provided a prescription to determine the r-function non-perturbatively along the 
integrable flow. 

A similar extension to OB turned out to be more subtle. Perturbatively, the computation 
is straightforward as it reduces to two copies of the c = 1 theory. These two copies can be 
thought of as either the left and right wave functions or, alternatively, the even and odd 
wave functions. 

Picking the even/odd basis, we hnd that, when R = {l/2)\/2a', the FZZT-brane deter¬ 
minant generates the appropriate Toda flow. Unfortunately, the flow is the same for the even 
and odd sectors. It is not clear if there is a better way to arrange the FZZT-brane so that 
^^Related discussion on integrable structure for minimal OB models can be found in [44]. 
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it only generates a flow in either the even or odd sectors. Perhaps a more serious problem is 
that even at i? = (1/2 )a/ 2^, we do not hnd that the resulting integral can be rewritten as 
a matrix integral. This issue arises from the fact that the integral is essentially two copies 
of the c = 1 integral, which required R = to be written in a matrix form. 

Non-perturbatively, which basis one chooses matters. In the case when the flux vanishes, 
the even/odd basis seems most natural as the problem continues to be manifestly factoriz- 
able. One can dehne two commuting Toda flows which generate separately even and odd 
deformations. 

However, when one turns on a flux, using the even/odd basis, one hnds that the OB parti¬ 
tion function is no longer T-dual to the partition function of the OA theory. If, however, one 
uses the left/right basis of [39], T-duality is restored. Unfortunately, in the left/right basis, it 
is no longer obvious how one should write down the Toda flow. Thus, to properly generalize 
the open/closed duality story to the OB case appears to require a better understanding of 
the integrable structure of OB. We hope to address this issue in the near future. 

There are a number of interesting open issues. One is to explore open/closed duality for 
Neumann FZZT branes. Unlike the macroscopic-loop operator for the Dirichlet case [32], the 
correct realization of the Neumann FZZT brane has not been found. There is an interesting 
proposal for representing FZZT branes using non-singlet sectors due to Gaiotto [53], but it 
does not appear to reproduce the correct disk and annulus amplitudes computed using world 
sheet techniques. 

An understanding of the Neumann FZZT branes might shed some light on our issues 
with the OB analysis since the T-dual of the Dirichlet FZZT branes in OB are the Neumann 
branes in OA. Unlike in the Dirichlet case, the Neumann brane acts as a source for winding 
modes. In the matrix model language, winding modes are encoded by the holonomy of the 
Wilson lines of the gauged quantum mechanics. The fact that there are two gauge groups in 
the OA model is related to the existence of the left and right sectors of the OB theory. Hence, 
if one can understand the integrable structure of the winding modes of OA in the presence 
of flux, it might suggest the correct integrable structure for the OB side. 

As a hnal point, we mention that the OA and OB models are simply two isolated points 
among a class of non-critical string theories, including affine type 0 theories and the type H 
theories. These theories are connected by a web of smooth deformations illustrated in [54]. 
It should be possible to explore the integrable flow and T-duality of each of these theories 
at the non-perturbative level, and to hnd a suitable generalization of the Imbimbo-Mukhi 
integral. 
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Appendix A: Conventions for Toda Integrable Structure 

In this appendix, we briefly review the machinery of the Toda hierarchy in order to fix 
our conventions. For a comprehensive introduction to the Toda integrable hierarchy, see 
[26,28,29]. 

Following the analysis of [23,24],we can the partition function of undeformed c = 1 model 
in terms of the Fredholm determinant. 


Z(/i) = det(l + e-^^A:_A'+), 


(A.l) 


where 



(A.2) 


Its log is therefore the sum over the eigenvalues of the operator, 

log(l + . 


(A.3) 


The orthonormal eigenstates, \E), can be expressed in the position basis as 



(A.4) 




(A.5) 


One then has 


log(Z(/i)) = {E\ log(l + = J dEp{E) log(l + (A.6) 


To compute p{E), one imposes a boundary condition. 


^^(A)=^^(A), 


(A.7) 
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which gives 


p{E) = ^ log A. (A.8) 

The log(A) term non-universal and can be dropped. By integrating by parts, one arrives at 
the standard expression, 

log(Z(^)) = - AI . (A.9) 

In the deformed theory, the K± operator is modihed; 

|if+/|(x_) = J f{x+), 

[A-_/|(x_) = /’[<ix_|e-“+'-/(x-) ■ (A.10) 

The eigenvectors of the operator K_K^ can be expressed as 



(A.ll) 

where 

(A.12) 

The operators VV± are called “dressing operators” and take the form. 


VVj. = W± = 

(A.13) 

where 

ch = . 

(A.14) 

They satisfy the relation, 

= 1 . 

(A.15) 

It is the relation (IA.15D. which, together with Theorem 1.5 of [26]. allows one to relate the 
computation of the deformed matrix model to the Toda flow. The structure of the Toda 
flow then allows one to compute 0(i?) and the v±nS as a function of ^±n’s, provided one 
specifies the initial conditions satished by the undeformed theory. Once the 0(£') is found, 
the partition function can be derived from (IA.9I1. 

Since we already know the wave functions which diagonalize K_K^ in 
theory, we can write the initial condition for >V± in the following form. 

the undeformed 

W±{Zn = 0,2_n = 0) = . 

(A.16) 


Here, the parameter a is introduced to allow for the ambiguity in the overall phase of the 
wave functions which does not affect any physical results. Indeed, different choices of a give 
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rise to the same flow [28]. We note this ambiguity here simply to warn the reader of the 
different conventions found in the references, with a = 0 and a = ±1 being the most popular. 

It is customary to parameterize W± as 

OO 

W^ = J2^f{E,z+,z.)u^^. (A.17) 

j=0 

Then, one can identify as being related to (p{E) in the parameterization of (jA.18jl ac¬ 
cording to 

w^{E, z+, z_) = e±td±a)0(i?) _ 

This form of W can be matched to the convention of [26] by identifying 

/OON 

Wj ° (s, z+, Z-) = Wj{E = i{s/R+ 1/2), z+, Z-) , 

OO 

0 \z+, Z-) = ^ diag[tt)j^ ° , 

j=0 

A = . (A. 19) 

With this choice, VV^ o) takes the form given in (1.2.9) and (1.2.10) of [26]. The sign and 
phase of z±n was chosen in (insD in order for the same z±n to appear on both side of this 
identification. 


Appendix B: Baker Akhiezer Function 


The Baker-Akhiezer function for the Toda lattice hierarchy is defined in terms of the dressing 
operators 


^ + s log AP"), 

Equivalently, one has 

^+(A) = T;(A) = . 

Defining the Lax operator, L±, and the Orlov-Shulman operators, M±, by 


(B.l) 


(B.2) 


(B.3) 
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M± = 

R 

it follows that the Baker-Akhiezer functions satisfy the differential relations, 

L^J,+(A) = A^UA), i-'tUA) = A®,-(A), 


(B.4) 


(B.5) 


d 


d 


M+fl-UA) = A—»yA), M_»;(A) = A—>1-;(A) 


d\ 


d\ 


(B.6) 


These Baker-Akhiezer functions are intimately related to the fermion wave functions of 
the double-scaled matrix model. Specihcally, 


V/TT 


(B.7) 


s=-ii(iE+l/2) 


Using the formal expression for the Baker-Akhiezer function, f)B.2jl . the orthogonality 
relation ()A.15jl . and the dictionary ()A.19fl . it follows that 


^*t(A+) j 


(B.8) 


The Fourier transform relation for the dual Baker-Akhiezer function (ESI, as well as the 
Bessel m and the sine/cosine transform for the c = 1 models, can be derived along 

similar lines. 

Appendix C: The Imbimbo-Mukhi integral for general n 

In this section, we show how the properties of the Baker-Akhiezer function determine the 
(n -|- 1) X (n -|- 1) Imbimbo-Mukhi integral in terms of the n x n integral. The computation 
is essentially the same as the 1x1 case, but the notation is a bit more complicated. It is 
useful to dehne 

n 

(c.i) 

i=l 


We then have the relation. 


T{-€n,Z-,S) 


(C.2) 


SO that 


r(-e(A -e„, 2 r_,s) = 
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Next, note that 



( Cri) ) -^+) 


(C.4) 



Substituting this into (nn gives 



X [ dXj6"'^+'^- - e(A!/^), s + 1) exp[z_(A_^^^) + (s + i?) log(Ai^^)], (C.5) 

J V 27r 

or, equivalently, using the notation, (j2.8njl . and shifting n —> n — 1, we hnd the recursion 
relation. 



(C.6) 


/ 


dM^ _y]^R 

_ p 71 n 


Zn-i{A,t- - e(M-^), s + i)e-’'^-iMr.)-is+i)iog(M^) 


X 


It is straightforward to check, by direct substitution, that solves (l(A6j) . 
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